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We study the edge and surface theories of topological insulators from the perspective of anomalies 
and identify a novel Z2— anomaly associated with charge conservation. The anomaly is manifested 
through a 2-point correlation function involving creation and annihilation operators on two decou- 
pled boundaries. Although charge conservation on each boundary requires this quantity to vanish, 
we find that it diverges. A corollary result is that under an insertion of a flux quantum the ground 
state evolves to an exactly orthogonal state independent of the rate at which the flux is inserted. 
The anomaly persists in the presence of disorder and imposes sharp restrictions on possible low 
energy theories. Being formulated in a many-body, field theoretical language, the anomaly allows 
to test the robustness of topological insulators to interactions in a concise way. 
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Topological insulators have attracted much attention 
in recent years due to their novel bulk and surface prop- 
erties. The bulk of these materials is insulating and 
characterized by topological indices which measure cer- 
tain twists in the band structure. The topological prop- 
erties of the bulk imply, via the bulk-edge correspon- 
dence, that the surfaces of these materials are necessarily 
metallic and support gapless excitations. The theories 
that emerge on the d — 1 dimensional surfaces of a d- 
dimensional topological insulator may be understood as 
" fractions" of the theories of " stand-alone" d — 1 dimen- 
sional systems [HH]. For example, the edge of an Inte- 
ger Quantum Hall (IQHE) at filling factor one, arguably 
the simplest and most striking topological state of mat- 
ter, may be thought of as half a spinless one-dimensional 
wire. 

The surface theories of topological phases have several 
important properties, which are clearly demonstrated for 
the case of the IQHE. Although IQHE edges border be- 
tween two insulating gapped phases, they appear to vi- 
olate a conservation law of the entire system, namely 
charge conservation. Indeed the low energy theory of 
each edge is invariant under the 17(1)— charge symme- 
try of the bulk. However, the application of an electric 
field parallel to the edges of a Hall bar leads to current 
flow across the bar. Consequently charges are exchanged 
between the edges and charge conservation, per-edge, is 
violated. This property is known as an anomaly and will 
occupy an important part of our discussion. In this case 
it is sometimes referred to as the Schwinger Anomaly [5] 
or a ID— Chiral anomaly [3]. 

While technically anomalies appear as a subtle cut- 
off scale effect, they are in fact present at all energy 
scales and have quite direct consequences [5]. For ex- 
ample the above Schwinger anomaly fixes the commuta- 
tion relations between the density operators on the edge 
[6] to a quantized non-zero value, while naively one may 
think that density operators at different points on the 
edge commute. Ignoring the anomaly in this case, one 
may arrive at the wrong impression that the edge con- 
ductance vanishes, as if the edge was insulating, instead 
of being e 2 /h. 



The identification of the anomaly associated with an 
edge or surface of a topological state is particularly easy 
for the case of IQHE, where charge conservation is being 
violated or for the spin Quantum Hall Effect (QSHE), 
where spin conservation is violated. Generalizations to 
other topological phases have been explored [2 [8] how- 
ever no anomaly was found to be associated with Z 2 
topological insulators (TI) in 2D or 3D. In this respect 
it is important to note an anomalous behavior, math- 
ematically similar to the SU(2) global gauge anomaly 
[3], associated with the 0(2N)/0(N) x 0(iV)-symmetry 
of the disorder-averaged action in the replica formalism 
[10] . However, since the symmetry is not a gauge symme- 
try, this anomalous behavior does not imply a symmetry- 
violation or an inconsistency in the theory. Thus, it re- 
mained unclear whether any symmetry becomes anoma- 
lous in Z-2 topological insulators. 

In this work we study topological insulators from the 
perspective of anomalies and identify a novel ^-anomaly 
that is associated with charge conservation on the bound- 
ary. The anomaly unifies the various aspects of topologi- 
cal insulators in a concise way through a field theoretical 
language and allows us to calculate new topological prop- 
erties of TIs. To formulate the problem we consider the 
quantum action of two decoupled boundary theories, cor- 
responding to two distinct edges or surfaces of a TI, and 
consider the cut-off scale as a bulk gap. We then dynam- 
ically insert a flux quantum so that an electric field along 
the boundaries is generated. Although the two bound- 
ary theories remain formally decoupled, we find that cer- 
tain quantities are exchanged via the anomaly. This is 
manifested in two ways. First, regardless of how adia- 
batically the flux insertion is carried out, the final state 
is always an excited states which is exactly orthogonal 
to the ground state. Second, even without any direct 
coupling between the edges, certain 2-point correlation 
functions involving creation and annihilation operators 
on different edges diverge rather than vanish. 

We begin by reformulating some well known results of 
the IQHE, at filling factor 1, in a way that allows a gen- 
eralization to TIs. Consider such an IQHE bar which is 
periodic in the x— direction and finite in the y— direction, 
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i.e., an annulus. The topological nature of the bulk is 
manifested by two chiral edge states with opposite ve- 
locity which appear on the two disconnected edges. The 
basic Hamiltonian describing each edge is simply 



Hedge = V (ihd x - eA x ) + /i, 



(1) 



where fi is the chemical potential and we allowed cou- 
pling to a gauge field (A x ). Below we study in detail 
the case where the edges are identical expect having op- 
posite velocities. Generalizations follow readily from the 
topological nature of our results. 

The spectrum of H e dge performs a spectral-flow as a 
function of flux (A x ) so that the states at A x = are 
exactly those at A x = h/(eL) displaced by a single state, 
where L is the edge circumference. Thus as one inserts a 
flux quantum adiabatically, exactly one electron is trans- 
ferred between the edges even though Hedge, on each 
edge, has a ?7(l)-charge symmetry. The spectral flow 
therefore implies the aforementioned charge anomaly of 
the IQHE edge. 

We find it convenient to introduce a formalism that 
treats both edges simultaneously, and express this 
anomaly as a 2D— chiral anomaly. To this end consider 
the action of both the edges following a multiplication 
of tjj by ia x 



s 



dxdT^ a [Sch}a,<T'i>o 



(2) 



S ch = (aihd T + iy)&x + v a a y (ihd x - eA x ) 



D 
Dt 



(3) 

where the a— spinor is associated with the edge index, 
ipjip are Grassman variables and a = l(i) for Euclidean 
(real-time) action. For simplicity we choose A x to be in- 
dependent of x. Furthermore it couples symmetrically to 
the edges since we do not allow time-dependent magnetic 
fields in the bulk. The action operator, S c h, acts in an ex- 
tended Hilbert space which includes the extra edge index 
and the time coordinate. Notably the U(l) symmetry as- 
sociated with charge difference between the edges is now 
reflected by the fact that {S c h,<7 z } = 0, which we from 
now on refer to as the chiral symmetry. 

Chiral symmetries are often anomalous [5] and the for- 
mer is no exception. Using Noether's procedure one finds 
that the chiral current is 



jch = ev ilja z aip, 



(4) 



where a — (a x ,a y ). A naive application of Noether's 
procedure will imply that the divergence of this current 
in space-time (d^jch,^) , with fj, = x,r, vanishes. This 
however is not necessarily the case in the presence of 
gauge-fields due to changes in the path integral measure 
[51 [11] . Instead one finds the following fundamental rela- 
tion of the chiral anomaly [5[ H2] , 



drdx(V jch) = Spectral flow 



drdxF (5) 



where — d t A x , A x is periodic, up to gauge transfor- 
mations, between — oo and oo (real-time) or —f3/2 to j3/2 
(Euclidean-time). By Spectral flow we denote the inte- 
ger number characterizing the displacement of the spec- 
trum of Hedge [A x (t)\ between t = — oo and r = oo. For 
example, Spectral flow = 1(— 1) implies that when fol- 
lowing an eigenvalue of ~H. e d ge [A x (T — — oo)] up to r = oo, 
one ends up with the upper (lower) consecutive eigen- 
value of the initial one. 

Equation ^ clearly reflects the physics of the IQHE: 
The charge transferred between the edges is related, in a 
quantized fashion, to the flux inserted into the annulus. 
This relation is easy to generalize for the spin quantum 
Hall effect, in which the chiral charge current Q is re- 
placed by a chiral spin current. However, for a generic 
2D TI there does not seem to be any local conservation 
law which is violated by an anomaly [7]. 

Chiral anomalies are also accompanied by the appear- 
ance of action zero-modes which do generalize to TIs and 
have important physical consequences. The zero modes 
of the current model can be understood by viewing S c h 
as a Dirac Hamiltonian on a 2D torus (spanned by (x, r)) 
with a uniform magnetic field. In this system a zero en- 
ergy Landau level appears with a degeneracy equal to 
the number of magnetic monopoles [13] . This relation 
turns out to generalize to any chiral operator coupled 
to a gauge field via the Atiyah-Singer (AS) and Atiyah- 
Patodi-Singer (APS) theorems [5] 



v = Spectral flow = 



drdxF, 



(6) 



where the analytic index v is given by v = DimKer[D] — 
DimKer[D^] and DimKer[D] (DimKer[D^\) is the di- 
mension of the zero mode space of D (D>). Note that a 
zero mode of D (D^) is also a zero mode of the action 
(Sch) with a z = 1( — 1) as can be verify explicitly from 
Eq. pi . Generalizations of all of our results to the // ^ 
case (non-hermitian action), are given in App. (AJ). 

Let us now generalize Eq. ([6]) to TIs. Consider a 2D 
TI cylinder with two distinct edges and again write the 
action for the two edges after performing the chirality 
transformation (ip — > ia x ip). This leads to a chiral action 
of the form 



Sch = (aihd T + i^Ox + <TyH e dge[A x ] = ( J)f Q 



D 



(7) 



where Hedge denotes the low energy Hamiltonian of a 
single TI edge, which we keep general. Provided that 
A x (t) — — A x (— t), the action is invariant under time 
reversal symmetry (TRS) 



TS ch T 1 = (iSyO- x P T )[S ch ] T (is y a x P T y 



(8) 



where s y = ^ ? ^ J is the Pauli matrix acting on the 
electron spin and P T flips the sign of r. Since T 2 = — 1 
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Kramer's theorem holds and each action eigenvalue is 
doubly degenerate. This again remains true also for the 
non-hermitian case in the sense that each Jordan block 
has an even dimension (see App. JA|)). 

The presence of Kramer's theorem for the action allows 
us to define a Z 2 analytical index similar to v. To see 
this let us first assume Hedge respects s z — symmetry and 
decouple the system into spin-up and spin-down compo- 
nents. When there is a single pair of zero modes, the two 
modes will have opposite s z and a z (see definition of T) . 
If we now break this spin symmetry, the zero modes and 
every other mode must remain degenerate due to TRS. 
The combination of the pairwise degeneracy of all modes 
along with the anti-symmetry of the spectrum around 
zero energy (implied by the chiral symmetry) allows the 
degeneracy of zero modes to change only in groups of 
four. This provides for the following topological index 

i>2 = DimKer[D] mod 2 = DimKer\D^] mod 2, (9) 

The TIs have an analog to the spectral flow of the 
IQHE edges in the form of " pairswitching" [2]. The 
"pair" refers here to the Kramer's pairs of edge states 
that characterize the spectrum at values of the flux that 
are time reversal symmetric, namely zero and one-half 
flux quanta. The term "switching" refers to the fact 
that states change their Kramer partners between zero 
flux and half a flux quantum. As we show below, this 
pairswitching behavior is related to the v%— index of the 
action via 



To relate Eq. ( 13 1 with Eq. ( 10 1 two observations 



v-2 = Pairswitching 



(10) 



The r.h.s. is equal to 1 (0) if the spectrum of H e dge[A x ] 
performs (does not perform) pairswitching as a function 
of A x (t). This generalization of the APS theorem relates 
the spectral motion characteristic of TIs boundaries with 
action zero modes in the presence electric fields. 

To prove the relation in Eq. |To| ), let us study the 
evolution of the action zero modes with r for A x (r) — 
If a zero mode, \(po(r, t), 0] t , with a z = +1 exists, 
the following equation must be satisfied 



iD(p = [ad T + H e dge[Ax(T)} - fl] tfo = 0. (11) 

with the boundary conditions, 

<po(r, r = -P/2) = -^(r, r = /3/2)e~ 2 ^ L . (12) 

To be concrete say we work with a — 1 (Euclidean 
action) and take an adiabatic limit so that (3 — > oo (the 
topological protection of the zero modes will later allow 
us to change these parameters). Provided that j3 is large 
enough the flux insertion can viewed as the adiabatic 
evolution of tpo in imaginary time. While it is not obvious 
that one may still use the adiabatic theorem in imaginary 
time, this turns out to be true under some limitations 
[15] . The (imaginary) adiabatic theorem then implies 

t)=J2 c n e- J "M^M-^te r), (13) 

n 

Tiedge[A x {T)]ip n (x\T) = [E n (r) -~ n]tp n (x; t). 



are important. First, pairswitching and TRS imply that 
there are an odd number of time points, denoted by r n , 
for which a state E m exists such that E m (T n ) = u and 
d T E m (r n ) is positive. Second, upon choosing c n = 5 nm , 
for each such point the adiabatic evolution yields 

ip {r, r) = e-^^-^V^r; r) + 0(e^ A ), (14) 



where A is the level spacing between two consecutive 
states of Hedge- In the adiabatic limit such states are 
localized near r„ and trivially satisfy the boundary con- 
ditions in Eq. (12 1. Therefore an odd number of solutions 
to Eq. ( 11 ) exist. Looking for zero-modes with a z = — 1, 
one obtains a variant of Eq. (11 1 with the sign of Hedge 
switched. Repeating the above analysis yields another 



odd set of zero modes. According to Eq. ^ this implies 
i/ 2 = 1 

In contrast, without pairswitching the edge spectrum 
(E m (r)) crosses with a positive slope, an even num- 
ber of times. Consequently, in the adiabatic limit, the 
e - / £ M- E ™( r )-M)_f ac t or m Eq. (13 1 guarantees an even 
number of zero modes for either D or yielding v 2 = 
according to Eq. 

Having established the stability of zero-modes and 
their relation to pairswitching, we wish to point at their 
physical consequence. To this end consider a TI in its 
ground state on an annulus, in which the flux thread- 
ing the hole is adiabatically varied from — $o/2 to $o/2 
($0 = h/e). We will now show that as long as the rate 
at which the flux is turned on is much smaller than the 
cut-off, the TI will end up being in a state that is exactly 
orthogonal to the ground state. Let us denote the ground 
state with — $o/2 by \gs) and consequently the ground 
state at $ /2 is G\gs), where G = exp l T J^+Wif-W. 
The state to which the system evolves after the flux is 
inserted is U\gs), with U being the time evolution op- 
erator. We examine the overlap (gs\G + U\gs) which, as 
shown in App. ([b]), can be written as the following path 
integral 



(gs\G^U\gs) = Z, 
Z = D[^V]« 



(15) 



dtdxiP[S ch {a(t))}iJ; + / dta(t)E gs , 



S 



where E gs is the many-body ground state energy, the 
limit of P ^> A -1 is assumed and we take a(t) = i for 
the period of the flux insertion (t G [—At, At]) and a = 1 
for t € [-fi/2,At], [At, P/2]. As shown in App. @, the 
boundary conditions are again those in Eq. (12). 



The path-integral expression can be formally evaluated 
[16] using the action eigenvalues to yield 



(gs\G^U\gs) = Det[S ch ] = II„/?A„ 



(16) 
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where A„ are the eigenvalues (or the Jordan eigenvalues) 
of (So?,,. Clearly in the presence of action zero modes this 
overlap vanishes. In the absence of zero modes and pro- 
vided that the flux insertion rate is sufficiently small, the 
adiabatic theorem implies that this overlap is one, up to 
phases. Thus the overlap serves as a sharp distinction 
between a topological and trivial insulator, formulated 
in a many body language. 

For the QSHE, the vanishing overlap follows from the 
change in spin on each of the edges, following the inser- 
tion of flux. For a generic TI, this behavior is much less 
obvious. Nonetheless in an extreme adiabatic limit, in 
which the flux insertion rate is smaller than A, it may 
be understood by elementary means. Consider the spec- 
trum of a TI edge and how it evolves with flux as it 
increases from — $o/2 to $o/2. At $o/2, all the states 
below the chemical potential are occupied. Following 
the pairswitching motion with flux, one finds that the 
ground state is degenerate at zero flux. This crossing 
of many-body states is not avoided due to TRS. There- 
fore, increasing the flux to $o/2 yields an excited state, 
orthogonal to the original ground state. 

Last let us examine the effect zero modes have on 
the chiral symmetry which, at least naively, implies 
the conservation of charge-difference between the edges. 
Consider the inter-edge two-point correlation function, 
Ghop = (J dxdtipaoSQip), where the average is taken with 
the chiral action. Since Ghop switches sign following a 
7T— chiral rotation, it should seemingly vanish. However, 
a direct calculation yields a different result. To show 
this, we include a source term S c h(m) — S c h — mipaoSoip, 
which physically amounts to coupling the charges on both 
edges, and find that 



G h op(m) = d m log 



Dk/^le - f dtdxt(:S ch (m)ip 

2v 2 



(17) 



Ghop appears as the time-averaged weak measurement of 
the inter-edge charge-tunneling operator. Physically this 
means that if one measures the charge-tunneling operator 
using a weakly coupled detector while selecting only the 
events in which the ground state came back to itself after 
the flux insertion- one will obtain a diverging value rather 
than a zero value. 

So far we have addressed the case of the 2D TI however 
generalization of this anomaly to weak and strong TIs 
readily follow. Weak TIs can be understood as layers of 
2D TIs and consequently, following its Z 2 — nature, the 
anomaly will either be absent or present depending on the 
number of layers. For strong TIs, pairswitching on the 
surface depends on two fluxes. Namely one flux induces 
pairswitching conditioned on the value of the other flux 
[18] . Consequently the anomaly, which follows directly 
from the pairswitching behavior, will be present or absent 
depending on whether this other flux is or $o/2. 

The Z 2 charge anomaly identified in this work, may 
facilitate our understanding of topological insulator in 
the presence of disorder and interactions. The anomaly, 
which persists also in the disordered case, guarantees 
a non-localized phase on the boundary of TIs. This is 
because anomalous theories cannot be massive or have 
only short range correlations at low energy [19 . The 
anomaly should also be present in the Non-linear-sigma- 
model (NLSM) descriptions of the disordered surfaces 
[TU1 [2171 [2Tj and may, through the anomaly matching con- 
cept, allow further investigation of these theories. Con- 
sidering interactions, it remains unclear whether topolog- 
ical insulators remain well defined. It is therefore inter- 



esting to see whether Eq. ( 17 1, which sharply distinguish 



= d rn log (n„/3[A„ + to]) = — + 0(m°, to 1 , ...) 



a topological surface from a trivial one, receives any per- 
turbative interaction corrections. 
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Appendix A: The anomaly in the non-Hermitian 



Interestingly, in the presence of zero modes, this two- 
point function diverges as the edges become more and 
more decoupled, rather than vanishes. For this rea- 
son we say that the chiral symmetry has an anomaly. 
Furthermore the appearance of 2v 2 as the coefficient of 
1 /to— pole, shows that this is indeed a Z 2 anomaly (see 
Eq. (|9|) so that coupling two anomalous systems yields 
a non-anomalous system. 

In quantum mechanical language, the anomalous cor- 
relation function can be expressed as a diverging weak 
value [T7]. To this end we note that following standard 
relations between coherent state path integrals and op- 
erators formalism one finds 



G 



hop 



At dt (Ggs{t)\[4 p^bouo 

-At 



h.c]\gs(t)) 



(Ggs(t)\gs(t)) 



(18) 



where \gs(t)) = Ul^gs), \Ggs(t)) = [U**\- l G\g8), U b a 
is the unitary evolution between t = a and t = b. Thus 



The chiral euclidean-action is an hermitian operator 
only for fi — and more generally, only when edge sym- 
metric terms are absent. Here we show that the results 
derived in the main text persist also when such terms are 
included. We first define the analytic index (y) for a non- 
hermitian action and prove its stability. Next a variant 
of Kramer's theorem is proven and used for showing the 
robustness of the Z 2 -analytic index (v 2 ). 

A general chiral action is given by 



S = 



L>_ 
D+ 



(Al) 



where D + and D_ are not necessarily hermitian conju- 
gates. The action, when viewed as a matrix, might be 
non-diagonalizable so that no spanning basis of eigenvec- 
tors exists, not even a non-orthogonal one. Nonetheless 
a weaker statement holds which is that any matrix can 
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be brought to a Jordan form (S) following a similarity 
transformation (S — PSP^ 1 ). In its Jordan form [22], 
the action becomes block diagonal so that each block ( J n ) 
has a single Jordan eigenvalue (A„) on its diagonal and ei- 
ther 1 or on the upper off-diagonal. All other entries of 
S are zero. Alternatively stated, a set of Jordan vectors 
exists (v n i) upon which the matrix (S — IX n ) is nilpotent 
so that (A — IX n ) Nn v n: i — 0, where N n < Dim[J„]. 

We extend the definition of v to the non-hermitian as 
followed 

v = Dim [J a {D_D+)\ - Dim [J (D + D_)} , (A2) 

where Jo (A) denotes the Jordan block associated with 
a zero Jordan eigenvalue of the matrix A. Let us 
show that for the hermitian case, the above definition 
is equivalent to v = DimK er[D^\ — DimKer[D], used 
in the main text. Note that for an hermitian operator 
Dim [j (DD^)] = DimKer[DD^} and similarly for the 
hermitian conjugate operators. Next note that if a vec- 
tor, v, is in Ker[D] (Ker[D^]) then it must also be in 
Ker[D*D] [Ker[DD^]). Furthermore, if v is in the ker- 
nel of Ker[DW] (Ker[DD^]) it must also be in Ker[D] 
(Ker[D^]) since D^Dv = implies v^D^Dv = \Dv\ 2 = 
and therefore Dv = 0. Consequently Ker[D] = 
Ker[D*D] and Ker[D*\ = Ker[DD^\ and the two defini- 
tions of v, for the hermitian case, are clearly equivalent. 

In the hermitian case each vector in Ker[D] (Ker[D'\) 
can also be used to build a zero mode of S with a z = 1 
(<7 Z = — 1) and thus the second definition of v relates 
better to the zero modes of S. For the non-hermitian 
case, the following weaker relation holds between v and 
zero modes of S 



Dim [J (S)] > 

To prove this last statement note that 

2 , I) I). 
1 D+D- 



(A3) 



(A4) 



and therefore Dim [Jo^ 2 )] > \v\ and using Jordan form 
we have that Dim [J (S 2 )] = Dim [J (S)}. 

Next we show that v, as defined above, is a stable 
number. Let v n ^ and A„ denote the Jordan vectors and 
eigenvalues of D + D-. Acting with D_ on v n ,i gives us 
Jordan vectors of D_D + , as one can directly show. Fur- 
thermore, D— generates a one to one map between the 
Jordan vectors of these two operators for which A„ ^ 0. 
Indeed, if there exists some linear combination v n , of v n> i 
such that D_v n = 0, then also D + D_v n = imply- 
ing A ra = 0. Consider now a small perturbation which 
increases Dim [Jo(D+D-)]. This implies that a Jordan 
block with some small A„ appears for and must 

be matched, through the mapping, to an equal size Jor- 
dan block appearing in D_D + . This later block can only 
appear from an equal change of Dim [Jq and 
consequently v is stable to any perturbation. 

Next we wish to show that the following variant of 
Kramer's theorem holds also in the non-hermitian case. 



Consider an action which obeys a Fermionic TRS sym- 
metry 



T SO = S T , 
O = -0 T , 
OO t = 1. 



(A5) 
(A6) 
(A7) 



where O in the usual case is is y . Then for such a matrix, 
S, each Jordan block has an even dimension. To prove 
this let us show how this symmetry acts on the Jordan 
form. Note that 



S 1 = [P^SP] 1 = P 1 S 1 [P 1 j- 1 = O 1 SO 



(A8) 



and by acting with [P 1 ] T and P T on the two sides of 
this last equality one finds 

S T = [P T ]- 1 T SOP T = [P t }- 1 O t P- 1 SPOP t . (A9) 

Denoting O = POP T one has that 

S T = O^SO, (A10) 



O 1 



-0. 



(All) 



Next we show that O is block diagonal on the basis, 
e n ,i, on which S and S T are block dia gonal . For each A ra 
block of size N of S, one can use Eq. ( A10 ) to show that 



o(5 J -\ n y = (s-\ n yo. 

Acting with e^ li from the left and using i(S— A„) 
one finds 



eiiPiS 1 -A„) iV W = 



(A12) 

N = 

(A13) 



for any vector w. Using the fact that S T — X n is invertible 
when constrained to the subspace of e m j with m =/= n, 
one may take w = [(S T — A„) JV ]~ 1 e m ^„ j and obtain 



e^iOe m ,j = Vm^ 



(A14) 



showing that O is indeed block diagonal on this ba- 
sis. Last using the fact that O is invertible and anti- 
symmetric we find that 



det[O n ] i= & Dim[O n ] e N e 



(A15) 



where O n is the n block of O. Therefore each block of 
the Jordan matrix must have an even number of states 
and a variant Kramer's theorem holds even in the non- 
hermitian case. 

The robustness of the ^2-index, for the non-hermitian 
case, follows readily from the definition 



u 2 = (Dim [J (5)] mod 4)/2. 



(A16) 



Indeed using the chiral and TRS symmetry of the ac- 
tion one can show that this quantity can only change by 
multiples of 4. 
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Appendix B: Overlap as a path integral 

Here we show that the partition function discussed in 
the main text is equal to an overlap of two distinct states: 
The ground state in the presence of a full flux quantum 
and the ground state to which we dynamically insert a 
flux quantum at a rate which is smaller than the cut-off 
scale or equivalently the bulk-gap. 

To establish this result we rewrite the overlap, up to 
phases, as 

(gs\G^U[gs) = Tr [PG^UP] = Tr [G^GPG^U] (Bf) 



P= lim e-Wl-^), 

0-H>O 

G — exp*T / d3KB V ,t ( a OV'(aO 



(B2) 

(B3) 
(B4) 



where H[A] is the full system Hamiltonian in second 
quantization coupled to a gauge field A, At is the dura- 
tion of the flux insertion process, E gs is a c-number equal 
to the ground state energy and therefore P projects on 
the ground state. The gauge field is given by A(t) = 
for t e [— At, At] corresponding to the insertion 



h 2LAt 

of a single flux quantum. 

Taking advantage of the fact that P is presented as a 
time evolution operator we may unite GPG^ , U and P 
into a single time evolution operator 



GPG^UP 



e - S%%dt a -\t)H[A(t)} ^ 



where for t e [-/3/2, -Ai] we take A = -e/(2hL) and 
oT 1 ^) = 1. For t e [At, /3/2] we take A = +e/(2hL) and 
a _1 (t) = f. In the remaining interval, t E [—At, At], we 

take A(t) = | ^TSt > a_1 (*) = *■ 

Using Grassmann calculus to take the the trace in Eq. 
(|B1|) we obtain 



(gs\G^U\gs) = / D[cc]e~ ec (-c|G t e _ J ^"/^ 



\c)=e-^\0), 



(B6) 
(B7) 



where |0) is the vacuum, a\ is a creation operator for some 
spanning basis of states indexed by i, cc — ^ZjCjCi, the 
limit (3 — > oo is assumed and summation over repeated 
indices is implicit. 

Next we follow the usual steps of constructing a path 
integral from Fermionic coherent states (see for exam- 
ple Ref. [23]). First the time evolution is partitioned 
into a product of N — > oo infinitesimal time evolutions, 
next Grassman resolutions of the identity in the form 
J"D[ip n 'iJj n ]e~^ n ^ n [ip n )(ip n \ are inserted between each of 
those, last using a n \ip n ) = tp n \ip n ) the operator H{a) ,a) 
is replaced by a Grassmanian functional. 

The only non-standard complication in this path inte- 
gral is the presence of the gauge transformation G. This 
turns out to be equivalent to a certain choice of boundary 



conditions as we now show. First we present the action 
of G on a grassman coherent state |c) through its action 
on the grassman variables 

G\c) = e- c * Ga ' Gt |0) = e- Ciff «°i|0), (B8) 
Ga\G^g %3 a\, 
Ga^G 1 = g^dj, 

so if i is the position basis, GajG -1 = e ^ 2Tra =i/ z - a | anc l 
9ij = Sije t27TXi / L . Following this we obtain 

G\c) = [eg), (B9) 
(cg\ = (c|G t = (0|e-°^ iff «. 

Boundary conditions for the path integral are deter- 
mined by the requirement that the discrete time deriva- 
tives generated by the time slicing process remain fi- 
nite. These time derivatives emerge from the —i^n^n 
factors coming from the resolution of the identity and 
the ■i/'nV'n+i coming from the overlap of adjacent coher- 
ent states. Focusing only on these elements, the path 
generates a series of the form 

e- Ec (-c 5 |^o)e-*' / ' (V'o|V'i)(-)e-^' / ' N (^|c> (BIO) 

_ e -cc-cg* tpo-ipoipo+tpoipl + i- ..)-i>NipN+i> JVC 

This series can regrouped to form a sum of derivatives 
either by grouping adjacent elements that share $ or el- 
ements that share ip (this ambiguity is a discrete version 
of integration by parts). Consequently we require both 
of the following expression to remain finite as we refine 
the time slicing 



-c(c+g'ip )-ipo('4'o~i('l) + {---)-i>N(ipN-c) 



(Bll) 



e (i>N-c)c+(-cg* -i(!o)-il>a + {4>o-4>i)-4>i + (---) +(ikN-i—ipN)i>N 

Identifying ip((3/2) with ipo and ip(—/3/2) with tpN the 
boundary conditions are 

M-P/2) = -5X/3/2) = -^{fi/2)g\ (B12) 

j>(-p/2) = -mmg T = -g$(P/2)- ( B1 3) 

which in position basis amounts to 

H-13/2, x) = -e- a ™' L i>{p/2, x), (B14) 

and its complex conjugate condition. The path integral 
expression for the overlap is therefore 



(gs\G^U\gs) = J B[^]e 



(B15) 



S = 4>d t ip - a- l (t)(H$,il>] ~ E gs ), 



with field configurations obeying Eq. ( B14 1 . Taking ip 
ipa(t), yields 



S = $a(t)dti/i-H$,rl>]+E ga , 



(B16) 
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In the main text we derived the zero modes associ- 
ated with the Euclidean chiral action (a(t) — 1), how- 
ever these also persist with the extra a(£)-factor. Indeed 
one can start from oto(t) — 1 gradually deform it, while 



respecting TRS as defined in Eq. to a(t). Such a 
continuous deformation cannot remove the topologically 
protected zero modes even though the operators associ- 
ated with S becomes non-Hermitian (see App. |A|). 
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